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ABSTRACT 


One of the main goals of the LAGEOS satellite mission is the 
detection of regional geotectonic movements. A parametric study with the 
intention to obtain the optimal baseline precision from dynamic solutions of 
laser ranging to LAGEOS is presented. The varied parameters are: length of 
reduced arc, number of tracking stations, data noise and rate, data biases, 
refraction errors, system efficiency, gravity model errors and errors in the 
value of CM. The baseline precisions are l-IO cm depending upon the set of 
parameters adopted. Geuerol principles obtained from the study are also 
presented. 
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1.0 INTRODUCTION 


Space techniques provide unique capabilities for measuring important 
parameters related to earth dynamics. They include laser ranging to satellites 
which is a precise tool tor the measurement of tectonic plate motions and 
crustal deformations over distances of hundreds to thousands of kilometers. The 
measurement and modelling of regional crustal deformations is an important 
subprogram of NASA's Crustal Dynamics Plan (NASA, 1979). The objectives of 
this program is to determine the physical mechanisms responsible for regional 
tectonic and geologic proccj^rs and to develop improved understanding of 
how regional phenomena fit to the fromework of global tectonics. 

The highly stable orbit of the LAGEOS satellite allows accurate 
estimation of the absolute geocentric positions of the afations of the loser 
tracking network (Smith et al., 1979a). The relative movement of sites, 
however, which is an important quantity for studies of regional crustal 
deformation, can also be determined through the estimation of relative station 
positions and intersite distances (baselines). 

Mobile laser systems with numerous site occupancies (the Mobile and/or 
Transportable Laser Systems) will play a very important role in the Regional 
Crustal Deformation Mi-osurement Program (NASA, 1979). The Crustal 
Dynamics Project Plan contains a detailed schedule for deploying the MOBLAS 
and TLRS facilities for regional studies, with western North America as the 
highest priority area. 

The plate velocities in areas selected for study in this program are of the 
order of a few centimr ters per year, so the regional-scale de format ional 
velocity measurements niust have an accuracy of a few millimeters per year. 
Such accuracies will permit comparisons of the regional motions with global 
ones, and will permit regional differences to be detected. Such accuracies are 
attainable with few-cent imeter accuracy in the measurement of intersite 
distances, within about ci decade of the beginning of the measurement program. 
For example, a baseline accuracy of 5cm from measurements made one year 
apart over a period ot 10 years, gives rise to a velocity uncertainty of 
0.5cm/yr. As the veloi ity uncertainties are directly proportional to the 

baseline accuracies, mu« h attention has been devoted to increasing the 
accuracy of the baseline determinations. 


The objectives of the present simulation study ore; 


a. fo obtain the general priruriples and improve our understanding of 
the dynamlcol reduction of satellite loser data for the determination 
of baselines of moderate length (200-500km), 

b. to identify the dominant sources of error ond recommend the 
optimal reduction method(s) that would minimize their effect on the 
baselines, and 

c. to obtain a qualitative assessment of baseline accuracies which can 
be expected g'ven our current and anticipated knowledge of the 
satellite's dynamic forces, data quality and geodetics of the 
reference system. 

The LAGEOS spocet raft was specifically designed for laser trocking aid 
is a completely passive satellite target. Orbiting at nearly 6000km altitude, 
LAGEOS is much less affected by the less well-known short wavelength 
features in the gravity field than lower geodetic loser satellites. Due to the 
spacecraft's spherical shape (60cm diameter) and its heavy weight (411 kg) 
LAGEOS is almost unaffected by the forces of the solar rodiation pressure, 
earth albedo ond air-drag. As a result, the orbit of LAGEOS moy be 
determined with high precision for time spans of several months. This 
characteristic of the LAGEOS's orbit plays o very importont role in the overoll 
improvement of the boseline accuracies. Typical Kepler elements ond rotes for 
LAGEOS tire given in Table I. 

Table I. Orbital Elements For LAGEOS 


Semi-major axis 12270 km 

Eccentricity .(XKi 

Inclination 109.84° 

Mode rate (C) 0.342°/day 

Argument of (lerigee rate (W) -0.24°/doy 

Period 225.5 min. 
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2.0 LASER DATA ANALYSIS TECHNIQUES 


The laser systems are highly accurate ranging instruments, which, 
while providing the rruist accurate means of tracking a near earth satellite, 
are not "all-weather systems." It is therefore quite difficult to schedule and 
subsequently achieve ideal tracking for numerous sites separated by great 
distances. It is of equal concern, in proctice, to be able to recognize the ideol 
tracking configurations, and likewise, find the optimol usage for the data 
which are actually obtained. While the latter is a far more common problem, 
the concept of ideol tiocking configurations worronts some discussion. The 
reduction of laser tracking of orbiting objects normally presents a dynamical 
modeling problem with numerous contributing forces (gravitational, rodiative, 
etc.), all of which are known to some level of imperfection. Consequently, the 
elimination of the reliance of force modeling altogether for the baseline 
solutions is a worthwhile gool. Such solutions, which rely solely on the 
geometry of the tracking and are taken simultaneously from the participating 
sites in sufficient number, provide unique relative positioning solutions for the 
respective sites which are no longer in any way dependent on accurate 
knowledge of the orbital position of the satellite. It is sufficient to know that 
all sites observe a common object simultaneously, irrespective of where the 
object is located. The data requirements, station configurations and other 
constraints of this geometric method are discussed in (Escobal et al., 1973). 
The geometric approach to the recovery of interstation distances has certain 
drawbacks. Firstly, tlie data requirements are severe and rarely, if ever, 
satisfied under the present deployment schedules for existing laser systems. 
Secondly, since the range data themselves exclusively define the solution to 
the station positioning, there is no external control (such as orbital dynamics) 
preventing the propagation of all measurement errors into the resulting 
solution. While in proctice, the laser systems perform quite well (Vonbun, 
1977), leoving this lotter concern of diminished importonce, it still cannot be 
disregarded if I to 3 cm baseline accuracies are desired. 

The objective of our simulations is the design and execution of on 
o ptimal solution using the available data to enhance the recoverability of the 
experimental objectives (precision baseline determination between the loser 


3 


network sites, or subst-is therein). Since, in all likelihood, dynamical methods 
will be tequired, a general discussion of this method within the context of 
precision station positioning is important. 

fven for L AGE OS, orbital determination at the 10 cm level requires 
extensive knowledge of several hundred terms in the geopotential, the 
atmospheric density, ballistic chorocteristics of the spacecraft, cross-sect ionol 
solar aspect area, lunar and solor masses and ephemerides, relativistic effects 
and others. Although satellite ephemeris error is not syrwnymous with station 
positioning error, as a general rule, ephemeris error during the tracking of the 
satellite is a significant error source for station positioning which uses the 
satellite’s orbit as a reference. 

Generally, over the years of experimentation with laser ranging, it has 
been deduced that by employing sliorter arc lengths to process subsets of the 
data, there is less propagation of the uncertainties of force model effects as 
errors in the calculated ephemeris of the spacecraft (Dunn et al., 1979). This 
con logically be extended to the geometrical case, where simoltcxieous sets of 
data poims ore processed individually. However, there are times when truly 
short arcs are inadequate for use due to poor data availability and the 
resulting inability to define the orbital reference system. 

Force model errors will be present, to vorying degrees, in all 
dynamical experiments. To minimize force model contomination, one can wait 
for improved force models, such as an improved harmonic coefficient model 
of the eorth's qeopoteutial, but such models are difficult to develop and may 
still prove to be of insufficient accuracy. Certoinly, the design of any 
experiment should have os an objective the minimization of the effect of 
dytiamicol errors in ttic meosured quantities. 

The San Andreas Fault Experiment (SAFE) is o cose in point 
(Sniith et al., 1979b). The baseline between San Diego and Mount Quincy has 
beer> successive. y measured numerous times within tlie past ten years. SiiK-i' 
the meosure of large scale plate motion was sought along the Son Andreas 


Fault, the experiment was successfully performed through on occeptonce of 
the presence of dynamical errors in e<x:h individual boseline measurement, but 
efforts were mode to moke this error constant for all subsequent bi>yearly 
recoveries. Therefore, groups of three successive passes of loser data on BE-C 
were used In each separate bi-yearly analysis. Geometry of the passes during 
the subsequent occupations was selected so that it was virtually identical to 
the originoi and the same force models were utilized. In this way, oil errors 
which are functions of geodetic position of the spacecraft manifest themselves 
in a very similar way for the 1972, 1974, 1976 and 1979 subset solutions. 
While the determined individuol baselines between Son Diego wid Quincy ore 
therefore biased, the relative fault nr)otion obtained through the dii erentlation 
of these results is a highly accurate meosurement, for the bios, year to year, 
has been mode nearly consvont. 
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3.0 


TLR5 BASELINl PRECISION FROM LAGEOS; A SIMULATION 


Tite particular focus of this error analysis is to assess the utility of the 
Transportable Laser Ranqing System (TLRS) for monitoring long term tectonic 
activity in the Western United States. The TLRS has now completed Its test 
and validation activities, and is octively deployed as part of NASA's Crustol 
Dynamics Program. It is a low energy laser operating using single 
photoelectrcn detection techniques (Silverberg, 1981). The system was built 
tor NASA by the University of Texas McDonald Observatory. The high mobility 
of the TLRS lies in its total containment within a single RV chassis allowing 
it to use crudely prepared sites with setup time requirements which seldom 
exceed a few hours. The analysis perfortned in this sfixJy is also applicable 
to other highly tnobile laser ranging systems, sixrh as the Compact Laser 
Ranging System (CLRS/ 1 LKS-IH which is undergoing development at Goddard 
Space Flight Center (Johnson, 1981). 

Our analysis has entered on estimating baseline accuracies for various 
TLRS deployment schemes In Caiitortiia. Figure I presents a map indicating 
the three TLRS sites wc> Irave considered. They are separated by distances of 
between JOO and TOO k'u. (Taselines in both North-South (in Figure I: the 

boseline from Site I to ‘nte ?) and Last-West (Site I to Site 3) directions were 
investigated. The MOItLAS systems, like t(x>se deployed at San Diego 
(SANDIF) and QUINCY have scheduled site occupancies of seveial months. 
They are therefore usucllv treated as fixed sites within our analysis. We have 
tried to address most of the practical problems related to TLRS site 
occupancies including length of stay requirements tor the TLRS, the 
requirements of the qlohnl laser network to support the TLRS activities, and 
various procedures tot minimizing the effect of dynamic errors in the 
reduction of the laser data tor baseline estimation. It is important to 
reiterate thot ttie deten linotion of intersite distances and not the geocentric 
station coordinates wert the experimental objectives. Consequently, baseline 
errors due to dynamic sources could be minimized although systematic errors 
in the stotion coordinnt! s were present. Station position errors were studied, 
fwwever, to provide som • insigtit into the characteristics of dynamical error 
propagation. 



In most of our simulalions we hove ossumed that the loser systems, be 
they TLBS or fixed instruments, have on efficiency of 50%. By this we mean 
that only half of all the visible Lageos possoges, complete and randomly 
selected, ore successfully acquired. This is in loir agreement with the history 
of system performances on the west coast of the United States. 

3.1 ORAN (ORBITAL /^ALYSIS) PROGRAM 

The simulotion study was occomplished through the use of o modified 
version of the ORAN program (Martin and Roy, 1972). This computer progrom 
simulotes a Bayesian least squares data reduction for orbital ond geodetic 
parometers. It does not process octuo! data. Through the generation of 
occurate normal equations, it has the capooility of computing the occuracy of 
the results if measurements of a given accurocy are available and processed 
in a least squares datu reduction program. ORAN is designed to consider o 
data reduction process m which a number of satellite data spens are reduced 
individually os well as simultaneously. 

The term arc refers to a specific data period over which a satellite's 
orbit is niegrnted; this dota solely defines the basis by which the satellite's 
position is adjusted. The effect of model uncertainties ond/or measurement 
errors are propagated into the set of basic porameters and into the estimated 
orbit over the arc interval. The effects of all error sources ore then 
statistically combined to produce a measure of the total resulting occuracy for 
both the orbit and ttw other adjusted parometers. ORAN was modified to 
accomn'odote the effrtt of a full voriance/covorionce matrix of o gravity 
field on the ailjrrsted parameters and to propogote the various errors into tfie 
baseline stotistics dircc tly. These modifications were necessary given the high 
correlation among certain coefficients of the low degree geopotential field as 
well as the need to thoroughly assess cancellation of errors into odjusted 
intersitr distances from stations which are in reosonable close proximity to 
one another. 

The error sources we have considered in our error analysis are 
summon/ed in Table ?. 



Table 2. Summary of TLRS Baseline Error Analysis Porameters 


Adjusted Parameter s; 

X, Y, Z, k, Y', i 

^|, A|, h| 

^ 2 " ^ 2 * ^^2 
(Implied Baseline) 


state vector of LACEOS at epoch 
for TLRS Site I 
for TLRS Site 2 


^Aodeled Parorreters: 


Dynamic Errors: 
Gravity Field 


100% of V/C motrix of GEM 9 
to 10 X 10 


GM 

Earth and Ocean Tides 
Ocean Loading 
Station Tidal Displocement 
Solor Rodiotion Pressure 
Measurement Errors ond Precision: 
System Efficiency 


I part in 10^ 

1% error In k2 
100% of effect 
10% error in h2 and /2 
0.5% of effect 


Measurement Intervol 
Ronge ^'loise 
Ronge Bias 

Tropospheric Refroction 
Elevation Cutoff Angle 
Reference System Errors: 
Fixed Station Positions 
Pole 4X, 4Y 


50% (complete passes randomly 
selected) 

I pt./sec. 

10 cm 

10 cm 

l%of effect 
200 

25 cm in eo< h coordirKJte 
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10 milliorc secoryfs 


3.2 


LAGEOS BASELINE SIMULATION RESULTS 


3.2. 1 G eneral Principles 

The Logeos baseline simulati<Kis were performed in cr> evolutionary 
process. Initially, certoin global properties of the baseline recovery problem 
were investigated. Subsequent simulations made use of the conclusions and 
principles of this eorlier global analysis permitting details of the TLRS 
baseline recovery to then be studied. 

The initial objectives of the study were to assess the minimum 
tracking requirements for the fixed laser stations of the global tracking 
network to support TLRS activities. A comparison was made between baseline 
recoveries for TLRS sites in which error sources were assumea to be constont. 
The tracking configuration was allowed to vary, however, so that: 

• In the first case, only the fixed west coast sites of Mount Quincy 
and San Diego were used, while 

• In the second case, six globally deployed sites (Figure 2) 
supported the orbit determination of LAGEOS. 

Table 3 summarizes the obtained results and shows that the estimated 
accuracy of the TLRS baselines was not dependent on global tracking support. 
Both cases yielded very similiar baseline error estimates. In this simulation a 
200 km (U-S) baseline was estimated between two TLRS sites within a 5 day 
arc length where all laser system were assumed to have 100% tracking 
efficiency. This level of efficiency, although unrealistic, was used so that the 
principle of global geometry requirements alone could be tested. Again, this 
simulotion showed that global tracking was not required. Theretore, for all 
subseq jent simulations we adopted a minimum requirement having only the 
west const fixed mobile laser sites supporting the TLRS. 
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FIGURE 2 











Table 3. (N-S) Baseline Accuracies (cm) 
Two TLRSs: 100% Efficiency (5 Day Arc) 



Case 1 

Case 2 

EF^ROR 

SOURCES 

5ANDIE 

& 

QUINCY 

FIXED 

6 GLOBAL 
STATIONS 
FIXED 

Positions 

0.1 

0.1 

Biases 

1.9 

2.0 

Refraction 

1.5 

1.5 

GM 

1.5 

1.4 

Gravity 

0.3 

0.3 

Noise 

0.1 

0.1 

Total (RSS) 

2.9 

2.9 
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The second conclusion from this initial simulation is summarized in 
Table k. The estimated accuracy (worst case) of the geodetic coordinates of 
the TLRSs are compared to their baseline accuracy estimate. The estimated 
accuracy of the determined baseline is approximately an order of magnitude 
better than the coordinates themselves. The dynamic errors, therefore, are 
shown to manifest themselves similarly in positional errors of nearby sites 
making highly accurate baseline estimations a realistic goal. This is similor 
to the conclusions of the SAFE experiment. 

We next sought on estimate of the effect of reduced tracking 
efficiency on the baseline error estimates. Table 5 presents the results of this 
phase of our analysis. Two problems were considered. Firstly, what was the 
impact of severe degradation of the tracking efficiency to 20% for the fixed 
sites along with the TLRS, (case l-b). In a second case, only the TLRS system 
was degraded (case l-c). It is opparent from these results thot the accuracy 
of the baseline is highly sensitive to the efficiency of the TLRSs and much 
less so to the efficiency of the fixed stations. The changing magnitude of the 
baseline errors due to gcopotcntia! uncertainties is responsible for this effect. 
However, as discussed later, this property of the problem has undergone 
extensive evaluation. 

Therefore, the general conclusions of this phase of our analysis are: 

o. global laser tracking is unnecessary for TLRS support if local 
tracking is available, 

b. the station coordinate errors are systematic, and while they are 
in the decimeter level, baseline precision is approximately an 
order of magnitude better, and 

c. when the TLRS efficiency is degroded below 100%, special 
consideration of geopotential errors is necessary. 
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Table k. TLRS Position and Boseline Accurocies (cm) 
Two TLRSs: 100% Efficiency, SANDIE AND QUINCY Fixed 

(5 Day Arc, Case I) 


EI^ROR 

SOURCES 

LATITUDI. 

LONGITUDE 

HEIGHT 

BASELINE 

Positions 

16.2 

17.8 

7.3 

0.1 

Biases 

6.4 

0.2 

15.4 

1.9 

Refraction 

2.9 

0.1 

5.9 

1.5 

GM 

27.6 

0.6 

41.8 

1.5 

Gravity 

7.4 

0.3 

10.0 

0.3 

Noise 

0.1 

0.1 

0.1 

0.1 

Total (RSS) 

33.6 

17.8 

46.6 

2.9 



Table 5. (N-S) Baseline Accuracies (cm) 

Two TLRSs: SANDIE, QUINCY Fixed (5 Day Arc) 
Variable Efficiency 


ERROR 

SOURCES 

Cose 1 

100% 

All 

Stations 

Positions 

0.1 

Biases 

1.9 

Refraction 

1.5 

CM 

1.5 

Gravity 

0.3 

Noise 

0.1 

Total (RSS) 

2.9 


Case l-b 

Cose l-c 

20% 

All 

Stotions 

20% TLRSs 
100% Fixed 

1,0 

0.6 

2.5 

2.4 

2.0 

1.8 

3.7 

2.1 

i6.7 

17.0 

0.2 

0.2 

17.4 

17.4 
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3.2.2 Optimisati on of Baseline Deler minotion for o Reolistic TLR5 Deploy- 
men4 i' nvironment 

Under the assumption of realistic tracking chorocteristics, the next 
phase of our analysis oddressed the following questions; 

a. whcil is the best reduction scheme for doto from a single TLRS 
having consecutive site occupancies? and 

b. what is the best scheme assumming two TLRSs ore providing 
laser ranging over the same time period? 

In this series t>l investigations we provided for a 50% tracking 
efficiency for all west roast sites porticipating in the trucking of LAGEOS. 
The length of stoy for the TLRS was varied, as was the redix;tion 
mettiodology. 

3.2.2. 1 One TLRS 


for the situation in which a single TLRS is occupying successive sites, 
nu(! tfn- determination ot the tKiseline t>etween these sites is the objective, 
numerous data reduction methods and varying lengths of stay for tlie TLRS are 
po:>sible. We hove attempted to evaluate the extremes of the suggested 
deployment schedule by varying the length of stay from five days (os in Bender 
ond Goad. 1979) to thirty days per site. Figure 3 is a pictorial synopsis of the 
configurations and data reduction methods which were evaluated. 

In the cases wh< re the length of the reduced arc was longer than 5 
days ,l|| and J]2 were odjusted in order to accommodate long period errors 
in the /orals of the geoiotential. This results in poorly determined individual 
coef fi( ients, but the res dting ephemeris will be relatively free of long period 
zonal fdasing. Since AGE OS has a high inclination we adjusted zonal 


(5 


TLRS DATA ANALYSIS CONFIGURATIONS 
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Ip 


base STATION (FIXED) 


coefficients of higher degree due to their insignificwt short period grovito* 
tionol perturbations which is not the case for the low degree terms. Short 
period effects will be somewhat altered nevertheless. For arcs of length 
shorter than 5 days, the adjustment of zonals is normally unnecessary due to 
the extensive ability of the orbital state vector adjustment to absorb these 
zonal errors. 

The data reduction schemes were evaluated for multiple short arcs 
versus arcs of longer length. Each configuration in Figure 3 is labeled (e.g., 
SB, 5S, etc.) for reference in the subsequent summary tables. Table 6 presents 
the breakdown of the contributing error sourc .*s in 200 km baselines (both N-S 
and E-W) for consecutive site occupancies of the TLRS. The baseline 
accuracies are computed by combining the largest errors from each of the 
error sources os obtained from four similar runs sp<»ming a 2 1/2 year period. 
They can be interpreted as approximately 2 cr accuracy estimates. 

Table 6 also presents a case in which a tracking efficiency of l(X)% 
has been utilized (solution 55/100%) for a five day arc length. This solution, 
consequently, has the maximum tracking geometry, the minimum errors from 
data noise, and is useful for comparison purposes. It provides insight into the 
relative degradation in baseline accuracy from sampling sources arising from 
poorer tracking efficiency. 

A discussion of the propagation of these error sources into the baseline 
accuracies is useful. It is important to understand that within a dynamical 
system, most error sources manage to corrupt the along track position of the 
calculated spacecraft ephemerides. However, the timing of the laser 
observations, and the minimization of the range observations itself at the 
point of closest (^prooch to the station, provides a meons for estimating an 
accurate period of the orbit notwithstanding the presence of these error 
sources. In other words, there is a strong dynamical control for the 
characteristics of the along track error (e.g., in the presence of tracking data. 
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they c«)not degrade secularly). When trocking both from the orbital ascent 
ar^ decent is available, along track errors tend to ultimately cancel. 
Systematic errors arising out-of-plone in the orbital ephemeris ore less subject 
to dynomicol control, and tend to tmve poorer properties for their cwKellation 
when estimating inter-station distances. Ttwrefore, longer arc lengths gen- 
erally are required for cancellation of cross track effects because the 
distribution of data longitudi '^ally must ochieve a balc^e. For the high 
inclinotion of the LAGEOS orbit, cross track is primarily in the longitudinal 
direction. The propagation of errors shown in Table 6 is discussed below. 

The fixed stati<»i position errors ore assumed to be uniformly 25 cm 
in e<K:h coordinate. However, one con deduce from the result shown in 
Toble 6, that the resulting along tr<x:k errors which ore in the direction of the 
N-S estimated baseline ore always smaller than their corresponding mcmi- 
festations for the E-W baseline, which is in the out-of-plone direction. 
Through an increase in the ore length, more error cmcellation becomes the 
cose, and the errors overall in both directions ore reduced. However, there 
is olways more tracking symmetry north to south about the TLRSs than eost 
to west. 

For data passes which ore not directly overhead, the best least squares 
accommodation of o ronge bias , if the pass is symmetrical from the northern 
to southern horizons, is an adjustment in the TLRS longitude. As a result, as 
shown in Table 6, adjusted E-W baselines ore more sensitive to bias effects in 
shorter arc lengths which hove a poorer data sampling in the longitudinal 
direction. However, for the thirty day arcs, good geometry is achieved and 
the N-S and E-W sensitivities become quite similar. 

A refraction error is similar to that of the range bias in that there is 
a non-zero mean error in range over the tracking interval. Like the range 
bias, this mean error can best be accommodated through an erroneous 
odjustment of the TLRS longitude. As seen in Table 6, for the shorter arc 
lengths, again, there is more sensitivity to this error source for baselines in 
the out-of-plone (E-W) direction. And for longer arc lengths, cancellation of 
this error source in both directions is achieved. 
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BASELINE ACCURACIES (CM) 


ONE TLRS: CONSECUTIVE SITE OCCUPATIONS 


ERROR 

SOURCES 

100% 

5S 

NORTH-SOUTH 

50% 

5S SB* 30S* 30/5 

30B* 

100% 

5S 

ss 

EAST-WEST 
50% 
5B* 30S* 

30/5 

3(»* 

POSITIONS 

0.2 

9.6 

1.6 

1.3 

IB 

U 

D 

12.9 

4.1 

3.3 

2.1 

3.3 

BIASES 

2.1 

3.9 

7.9 

2.0 

2.5 

2.0 

0 

8.0 

3.5 

1.8 

1.3 

1.2 

REFRACTION 

1.7 

2.2 

2.3 

1.6 

1.9 

1.7 

0.5 

3B 

IB 

0.7 

0.9 

0.6 

GM 

1.7 

3J 

5.2 

1.4 

2B 

1.1 

SI 

18.7 

14.5 

5.0 

1.8 

1.1 

GRAVITY 

1.3 

14.1 

10.3 

5.3 

5.4 

4.1 

u 

13.3 

6.5 

3.7 

4.1 

2.1 

NOISE 

0.1 

0.2 

0.2 

0.1 

0.1 

0.1 

0.1 

OB 

0.2 

0.1 

0.1 

0.1 

TOTAL! RSS) 

B 

18.0 

12.2 

6.2 

B 

5.1 

2.9 

27.9 

16.9 

7.3 

5.2 

4.3 


S77 

S76 

C76 

C?6 

S76 

C76 

C33 

S77 

C76 

C77 

C77 

C76 


C77 

C76 

S22 

S76 

S77 

S76 

C76 

C76 

S76 

S77 

S76 

S76 

GRAVITY 

S97 

C77 

S76 

C66 

C78 

C44 

C22 

S76 

S22 

S76 

C30 

C77 

FIELD 

CB7 

C22 

C77 

S22 

Gs6 

S77 

C66 

C66 

C3.3 

S44 

C33 

S22 

TERMS 

C20 

C44 

C66 

C77 

See 

S22 

S76 

C30 

C66 

C22 

S22 

S41 

IN ORDER 

C40 

S22 

C44 

See 

S97 

S44 

S22 

S97 

S21 

S31 

C76 

See 

OF 

$66 

S77 

C21 

S77 

C77 

See 

C20 

C20 

S31 

C76 

C66 

S77 

IMPORTANCE 

C76 

S44 

S21 

S86 

S22 

C66 

C40 

C40 

C21 

C43 

C22 

C21 


PB5 

C97 

C41 

C44 

C22 

C77 

S21 

S31 

See 

Gee 

C44 

S31 


C22 

S43 

C97 

C33 

QB6 

C21 

See 

S21 

C77 

S97 

S77 

C22 


C66 

S97 

C33 

C85 

S43 

Ceo 

C53 

C77 

S77 

C31 

C97 

C97 


S76 

C43 

S3I 

CB7 

S44 

C20 

C30 

S41 

S41 

S33 

See 

S21 


*Jl ^ . J)2 adjutttd to account for taculw affects. 


TABLE 6 
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The situation for errors in GM is more complicated. As alluded to 
earlier, the orbital period (mean motion) con be inferred directly from the 
loser system timing, the point of closest opprooch wkI the adjustment of the 
semi-mojor axis of the satellite. GM is o*sumed to be held constwit in our 
simulations, and ce o consequence of Kepler's third law, the fractional error 
in GM is approximately equal to one third of the error in the odjusted semi- 
major axis. This in tinm, scales the size of orbit to the adopted (erroneous) 
value of GM. Therefore, since the rcrtge data hos not undergone a similw 
scoling, the oppwent errors in range ore again similar to those of the roige 
bias. As seen in T^le 6, the s<xne principles os those opplic<^le to range 
biases apply for the propagation of the errors due to GM. Systematic 
longitudinal errors result and require longer arc lengths for cancellation. 

The dynamical error iources arising from the uncertainty of the 
geopotential field dominate the degrodation of the estimated boselines in 
Table 6. A bredcdown of the individuol sf^rical harmonic coefficients within 
the field reveals that orbital resonance with ms6 terms are the largest single 
error source (see Appendix A). However, terms hoving large m-daily effects 
are also significwit (see Apperxlix A). The odjustment of the zonal harmonics 
hos greatly reduced the Impact of zornil hormonic uncertainty for the longer 
arcs. The gravity model error results in a complex degradation of the along 
track ep^merides accuracy of LAGEOS which is three to five times worse 
thon the corresponding cross trock errors. However, these errors ore not 
symmetrical in either north to south or east to west directions. Given the 
magnitude of the along track vs. cross track effects, gravity hos more severe 
consequence for N-S bosellne odjustnrwits th<r» those E-W. However, with 
improved sampling through o lengthening of the ore, cancellation of the 
gravity error is apparent, but residual error is still quite substantial. 

The influence of doto noise is directly proportional to the square root 
of the number of observations. Therefore, long arc lengths, with more data, 
show improved noise-only baseline accuracy estimates. 
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T<rt>le 7. Baseline Error (cm) 


One TLRSt Consecutive Site Occupations; 
50% Efficiency (70 Day Arc) 


ERROR SOURCE 

ERROR MAGNITUDE 

BASELINE ERROR 

Solar Rodiotion 

0.5% in Cp 

(15.00)* 

Earth Tide 

1.0% in l<2 

0.5 

Geometric Tide 

10% in h2 and ^2 

0.2 

Pole 

0?0l in X ond Y 

1.4 

Oceanic Loading 

100% of Effect 

0.0 


♦Requires Adjustment 
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It is (^>parent from this phose of our analysis that lon 9 «r ore lengths 
yield improved baseline accuracies for a single TLRS having multiple 
consecutive site occt^Mncies. Therefore, based solely on the <^ve onolysis, 
we recommend TLRS site occupancies of at least thirty doys if a sirigle system 
is employed and boselines between TLRS locations ore the experimental 
objectives. 

The other error sources (rx>t menti<xied in Table 6) need to be consid- 
ered also for tlwse longer arc lengths. Table 7 presents this information os 
the worst error, in either the E-W or N-S directions, arising from these os yet 
unmentioned error sources. As is evident from this table, a solar radiation 
pressure coefficient needs to be odjusted for ores of 70 days length. However, 
all other error sources are relatively insi^ificont. From further simulations 
we have found t»«t on adjustment of a solar radiation pressure coefficient has 
no influence on the baseline errors arising from any other source. 

3.2.2.2 Two TLRSs 

Our simulations have also deoit with o tracking configuration where 
there ore two TLRSs avoiloble. If both systems are deployed ot the some time 
and scheduled to track LACEOS during the seme working hours, then a new 
data reduction methodology can be attempted. One con limit onolysis of the 
TLRS data sets to those observations which ore simultaneously avail<rt>le from 
both sites. This is grophicolly displayed within Figure U where the 
simultaneous data is token over the interval defined between points 2 and 3. 

This definition of simultaneity is strict, and by this we meon that: 

a. observations must be temporarily matched to within reosonoble 
limits, ond 
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2k 


FIGURE 


b. the simultaneous passes in the arc must contain approximately 
the same number of observations (allowance can be made for a 
small number of segments of passes to have fewer points). 

The data from the fixed stations are not subject to these restrictions and all 
of their available data is used in all cases. 

Our interest in this q>prooch lies in our belief that the fixed stations 
can sufficiently define the orbit in the vicinity of the TLRSs. Further, by 
limiting the TLRS data set to times of strict simultmieity, we ensure 
maximum cancellation of the errors in the estimated 200 km TLRS baseline. 

In this series of simulations we have again assumed a tracking 
efficiency from all systems of 50%. Consequently, one would expect strict 
simultaneity to occur some fraction of this time. We have further degraded 
the efficiency of the systems by assuming that only 50% of the observations 
within the pass itself are successfully obtained. Therefore, the passes are now 
incomplete, and simultaneity needed to be assessed on a point by point basis. 

Obviously, if this approach was shown to be optimal, long arc lengths 
are not required to provide cancellation o' error source. This is the result 
indicated through a simulation of a five day arc length. 

Table 8 intercompares the baseline accuracy estimates obtained from 
a configuration having tracking from two TLR5s for five day arc lengths. In 
cases 2-a|s| and 2-a£ (the N denoting a N-5 baseline, E for E-W), a tracking 
efficiency of 100% was used. As before, this test case was used as a standard 
by which degradation (or in our case, possible improvement) could be assessed. 
Cases 2b|vj and 2b£ used all data provided by the algorithm of 50% tracking 
efficiency discussed earlier where both passes and observations are removed. 
Coses 2 c|m and 2c£ use subsets of the data from 2b, but for the TLRS only 
those observations adhering to the simultaneity requirements are used. 
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BASELINE ACCURACIES (CM) 

TWO TLRS*S: SIMULTANEOUS DATA CONFIGURATIONS (5 DAY ARC) 




NORTH-SOUTH 


EAST- 

-WEST 



1 50% of poinbi; Incomplete penes 


50% of points; Incomplete pMses 

ERROR 

All 


Ban Statiwu 

All 


Ban Stations 


Pot- 

All 

All Pastes 

Pot- 

All 

All Passes 

SOURCES 

sibit 

Penes 

TLRS'sSimult. 

sible 

Pastes 

TLRS's Simult. 


PanM 


wily (9 PMtes) 

Passes 


only (12 Passes) 


2 »n 

a»N 
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2 bg 

2 Cg 

POSITIONS 

0.1 

0.6 

0.1 


0.3 

0.1 

BIASES 

1.9 

1J 

0.2 


0.9 

0.3 

REFRACTION 

1.5 

1.5 

0.4 

0.1 

0.6 

0.5 

GM 

1.5 

1.5 

0J5 

2.1 

5.0 

0.7 

GRAVITY 


10.9 

OB 

■ 

3.6 

0.3 

NOISE 

0.1 

0.2 

0.3 

■ 

0.1 

0.2 

TOTAL (RSS) 

2.9 

11.3 

1.1 
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1.0 



*76 

S6 

^33 

S2 

S2 

GRAVITY 

S 77 

<=76 

s? 

Se 

*76 

Sb 

FIELD 



s„ 

c»« 




32 

44 

77 

66 

66 

20 

TERMS 

^22 

*77 

^66 

S2 

S 7 

So 

IN ORDER 

S? 

^44 

S 7 

S2 

^33 

CM 

OF 

S31 

CM 

S2 

S3 

S3 

*21 

IMPORTANCE 

®32 

^5 

So 

So 

Sb 

*41 


Se 

Si 

<=<0 

*31 

<=43 

*76 


TABLE 8 
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One can see that in all cases, a five day arc containing two TLRS sites 
yields baseline accuracies which are superior to those obtained from two 
successive five day solutions (solution 5S in Table 6) from a single TLRS 
occupying consecutive sites, t^wever, when the data set of cases 2b are 
employed in their entirety, the cancellation of errors in the baseline is not 
very successful. This is primarily due to the along track errors (again 
principolly resonance) arising from the uncertainty in the geopotential. 

The results for the solutions employing strictly simultaneous data sets 
(cases 2 cn and 2c£) are extremely encouraging. The expected cancellation of 
unmodeled errors in the determination of the baseline has occurred to a very 
high level. Gravity errors are now reduced to a single centimeter error in 
baseline. Although orbital resoncxice error is still dominant, its effect on the 
determined baselines is now reduced by an order of magnitude. 

Table 9 presents the estimated effects of all other considered 
parameters in the simultaneous data solutions: 2cfsi 2c£. Obviously, none 
of these remaining errors are significant beyond the I cm accuracy level. 
Again, Table 9 presents the worst incidence of error, in either the N-S or E-W 
baseline. We therefore conclude that based upon these results, the 
simultaneity method is optimal for analysis of data from two TLRSs. Baseline 
accuracies of I to 2 cm are possible through the implementation of this data 
reduction technique for those baselines of 200-500 km length. 


3. 2.2.3 One TLRS: Application of the Simultaneity Principle 

The principle of employing strictly simultaneous data between adjusted 
stations has some application for a situation in which only a single TLRS is 
available. However, to do so, one of the previously fixed MOBLAS sites must 
now be treated as though it was the second TLRS. MOBLAS must be adjusted 
and only its data which is strictly simultaneous with that of the TlRS are 
used. The baseline of interest in this case is now the line between the 
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Table 9. Baseline Error (cm) 


Two TLRSs: Simultaneous Data Configurations 
(5 Doy Arc) 


ERROR SOURCE 

ERROR MAGNITUDE 

BASELINE ERROR 

Solar Radiation 

0.5% in Cr 

0.0 

Earth Tide 

1.0% in l<2 

0.0 

Geometric Tide 

10% in h2 and ^2 

0.2 

Pole 

OVOI in X,Y 

0.2 

Oceanic Loading 

100% of Effect 

0.0 
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MOBLAS and the TLRS, and not those lines between successive TLRS site 
occupancies. To minimize the errors between TLRS baselines exclusively, ttw 
methodology detailed in Section 3.2.2. 1 is still superior. 

Accurate baseline determination between the TLRS and the MOBLAS 
is possible through the employment of the simultaneity principle. This Is 
presented in Table 10, where SANDIE (Son Diego) is now treated as on 
adjusting station along with the TLRS at site 3 (from Figure I). The baseline 
we are now discussing is 413 km in length, over twice the length of <my 
previously considered line. Again, as in coses 2c (in Section 3.2.2.2) the 
tracking efficiency level employed for all stations is 50% for both passes cvid 
individual points within the pass. There is now only a single fixed station, 
QUINCY, for the first case considered, which is case 3a on Table 10. All data 
from QUINCY allowed by the algorithm are used while SANDIE aid the TLRS 
supply strictly simultaneous data. 

In case 3a the orientation of the baseline in three dimensioial space 
is inferior to the previous situations using two fixed stations. This distortion 
increases the error propagation into the baseline from cross track effects such 
as those arising from range biases and refraction errors. However, these error 
sources have been unduly pessimistically modeled throughout our simulations, 
and are subject to improvement. What is more significant, however, is the 
relative insensitivity of this approoch to gravity error sources even though 
resonance error dominates. 

In case 3b, a second fixed site located in Mexico, is introduced. The 
orbital plane again becomes more stabilized, and the cross track error sources 
diminish in the baseline statistics. We therefore conclude that the 
simultaneity principle within the framework of a dynamical data analysis 
warrants consideration if the baselines between a MOBLAS and the TLRS sites 
are in the direction of tectonic interest. If this is so, baselines of 2-3 cm 
accuracy between the TLRS and the MOBLAS sites are achievable. 
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BASELINE ACCURACIES (CM| 

ONE TLRS: SIMULTANEOUS DATA CONFIGURATIONS (6 DAY ARC) 


ERROR 

SOURCES 

SANDII 

(413 

50% of poini 
Oulncy All Pmw 
SANDIE&TLRS 
Simultmoout Only 
(9 tawt) 

3a 

SITE NO. 3 
Km) 

i; Incompitt* pantt 

Ckiincy & Maxico 
AIIPmim,SANDIE 
& TLRS SimultwMout 
Only (9 Pmw) 

3b 

POSITIONS 

0^ 

0.2 

BIASES 

3.7 

1.7 

REFRACTION 

3.1 

1.4 

GM 

0.9 

1.1 

GRAVITY 

1.0 

0.7 

NOISE 

0.3 

0.2 

TOTAL (RSS) 

5.0 

2.6 


*77 

Se 

GRAVITY 


^66 

FIELD 

Se 

*76 

TERMS 

Se 

*77 

IN ORDER 


®66 

OF 

*97 

*44 

IMPORTANCE 

^21 

^32 

j 

^22 

Se 


TABLE 10 
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A.0 CONCLUSIONS 


In the present study we elaborated on the expected TLRS baseline 
accuracies in the Western United States using the LAGEOS satellite. The 
conclusions we could draw from this error analysis regarding the accuracy of 
baselines of moderate length (e.g., 200-500 km) are given below. 

a. Global laser tracking is unnecessary for TLRS support if local 
tracking is available. 

b. The station coordinate errors are systematic, and while they are 
in the decimeter level, baseline precision is approximately on 
order of magnitude better. 

c. When the TLRS efficiency is degraded below 100%, special 
consideration of geopotential errors is necessary. 

d. For baselines between consecutive TLRS site occupations; 

• 30 days on site is necessary as baseline accuracy is 
directly dependent upon efficiency of TLRS, 

• Above requirement may be lessened if the gravity field is 
improved, and 

• Long arc reductions slightly favored over multiarc reduc- 
tions in which case the 2cr accuracy of the baseline is 
better than 5 cm. 
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e. For baselines between bose station and TLRSj 

• If both adjusted cvid only simultaneous passes allowed in 
the solution the 2e accuracy of the baseline is better then 
I part in 2*10^. The reduced arc can be as short as 5 
days and the number of strictly simultaneous passes could 
be 6-8 {not necessarily complete). 

It must be kept in mind that some of the above conclusions could possibly be 
radically revised after the first couple of years of observing, as a result of 
improved knowledge of previously poorly known quantities. 
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APPENDIX A 


Resononc e 

Resonance is a common effect experienced by near-earth orbiting 
objects. Basically, resonance is a short period longitudinal-dependent term of 
the geopotentiai which monifests itself os on excessively large long period 
perturbation on the orbit due to the commensurability of the satellite's 
motion with the earth's rotation. 

For most orbital arc lengths which exceed even a few revolution 
lengths, reson<ax:t presents itself as a signific(^t problem. However, the 
terms in the field which resonate with the orbit generolly have a single 
significant contribution to the orbitol motion. LAGEOS, which has a me<n 
motion of nearly six revolutions per doy, is resonant with m=6 terms of the 
geopotential, C,S(7,6) being domincnt. The contribution of all m=6 terms, 
however, is almost entirely due to these resonance effects. Table A shows 
all significant perturbations arising from C,S(7,6) on the LAGEOS orbit. It 
is evident that the resononce effect arising from these terms is several 
orders of magnitude larger than any resulting short period or m-daily 
perturbation. Errors resulting from imperfect knowledge of the resonance 
harmonics con be arcommodated through the odjustment of a select set of 
resontBice terms. While the individual coefficients obtained th:ough this 
adjustment moy suffer from severe aliasing, the resononce error is minimized 
in the orbitol trojectory. 

M-Doily ond Short Period Terms nf the Geopotential 

The dominant source of low degtee and order information (/< 8) in the 
geopotential comes from the m-doily perturbations which manifest them- 
selves on all of the orbits used in grovitationol model recovery. The m-doily 
terms orise from short period perturbations which are independent of (i.e.. 
averoged over) the mecn anomaly. However, the accommodation of errors 
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TABLE A 

ESTIMATED ALOMG TRACK FERTURBATIOMS 
FROM TRE DOMINAMT RESONANCE TERMS ON THE LACEOS ORBIT 


HARMONIC 

CONSTITUENT 

FREQUENa 
(eye /day) 

ESTIMATED 
ALONC TRACK 
PERTURBATION 

In) 

GRAVITATIONAL 

FAMILY 

i E E a 

7 6 3 0 

0.37 (2.66 day 
pariod) 

28.28 

ahallow resonance 

7 6 2 0 

13.15 

0.13 

short period 

7 6 2 1 

19.53 

0.06 

short period 

7 6 4-1 

18.78 

0.08 

short period 

7 6 4 0 

12.40 

0.26 

short period 

7 6 3 1 

6.76 

0.30 

short period 

7 6 5 0 

25.17 

0.06 

short period 

7 6 6 -1 

44.33 

0.10 

short period 

7 6 6 0 

37.94 

0.12 

short period 

7 6 7 0 

50.71 

0.06 

short period 

7 6 3 -1 

6.01 

0.18 

a - dally 

7 6 4 1 

6.01 

0.18 

a - dally 


Based on 

Normalized value: ^^^2 

■>7.6 • * *7.6' ) ■ 


*From CEM9 (Larch ct al, 1977} 
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from imperfect knowledge of the m-doily effects is mode very complex due 
to the presence of a rich spectrum of other short period perturbations arising 
from the some low degree and order harmonic coefficients. As an 
illustration, Table B presents all of the significcnt orbital perturbations on 
LAGE05 produced by C,S(4,4) and C,S(6,4). The amplitude of each 
perturbative frequency consists of a linear combinatian of the effects of 
terms of the some order (m) having consistently either odd or even degree 
as shown in Table B. All of the frequencies of the C,S(4,4) terms ore found 
in terms of C,S(6,4) although the reverse is not true. These very same 
frequencies will r'oo be found in every higher even degree, 4th order term. 
However, the re itive importance of each frequerKry (e.g., percent of 
dominant constituent in Table B) is a function of the orbital elements of 
the satellite and also a function of the harmonic's degree for a given order. 
Although C,S(4,4) has only two harmonic constituents producing nearly its 
total perturbation ( / mpq=(4,4,2,0) and (4,4,2,!)) on LAGEOS, C,S(6,4) finds 
these two frequencies having much less importance. As shown, C,S(6,4) 
manifests itself on LAGEOS in five constituents which ore at least 40% of 
its dominant effect. With the eccentricity being small, terms of odd degree 
will have frequencies very similar to those arising from the even degree. 
This further complicates the gravity modeling problem. The same type of 
behavior is exhibited by the entire low degree and order field. 

As a consequence of these numerous frequencies and the large number 
of the resulting linear combinations of the harmonic coefficients a simple 
adjustment of a select set of coefficients presents a difficult problem. This 
is especially true given the nature of laser tracking, with the data acquisition 
problems alluded to in Section 2.0. Laser systems are not all-weather 
instruments, and consequently, the data which is obtained can vary 
dramatically from site to site, and even day to day for a given site. 
Therefore, there are significant problems which arise from incomplete 
sampling when one is attempting to uncouple a large number of mismodeled 
short period orbital effects from station coordinate errors if the geopotential 
coefficients are adjusted. 
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TABLE B 

SIGNIFICANT GRAVITATIONAL PERTURBATIONS FROM C,S(A,4) 
AND C,S(6,4) ON THE LAGEOS ORBIT 
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A phenomenon similar to resonance can arise between the tracking 
from 0 given station and certain frequencies of geopotential error. For 
example, the LAGEOS arbit is visible to a given mid-latitude site at 
approximately 12 hour intervals (i.e., one sees the satellite's ascent, the 
earth rotates beneath the orbital plane, and 12 hours later, the satellite's 
descent is viewed). In our LAGEOS simulations one finds a large latitude 
error resulting from m=l geopotential errors. Since m=l terms have m-daily 
effects which have a one cycle/day frequency, the net effect is the along 
track position is leading at a given time to be followed by a corresponding 
log 12 hours later due to the errors in these m-daily m= I terms. Considering 
the high inclination of LAGEOS, along track effects are nearly totally in the 
latitudinal direction. This orbital error therefore, can be directly absorbed 
(i.e., the range errors can be minimized) by moving the station's latitude an 
appropriate amount. When the data sampling is incomplete, the mani- 
festations of the mismodeled low degree short period effects, becomes 
largely a function of the sampling itself, and becomes unique to each specific 
arc. This problem is greatly magnified and the interpretation of simulated 
results becomes far less uncertain, when the low degree part of the 
geopotential is allowed to adjust. The quality of the geopotential adjustment 
and the net resulting baselines can vary significantly if differing random 
tracking configurations are employed. The determination of a complete 
(albeit truncated) geopotential requires tracking from numerous satellites 
taken over long periods of time. Only through the analysis of this type of 
data set can one uncouple the linear combination of gravitational effects into 
well determined individual harmonic coefficients.* Data spanning a full 
apsidal period is normally required to separate odd from even degree effects. 


*Tailored fields have the correct linear combinations to produce the correct 
amplitudes and phases on a given orbit but not necessarily the correct values 
for the individual gravitational coefficients. 
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